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Lattice simulations of B ^ ;r/v decays are problematic in the low region when the pion has to 
have large momentum in the lattice frame for a B meson at rest or moving slowly. The moving- 
NRQCD formalism provides a way around this by giving an accurate discretization of a b quark 
moving at arbitrary velocity v, and therefore a B meson moving at high momentum in the lattice 
frame. Here we show results from simulations of the moving-NRQCD action complete through 
0{l/M) coupled to the asqtad action on a MILC fine ensemble. We show B, B* masses as a 
function of v and demonstrate how to determine the meson kinetic mass and the renormalization 
of V non-perturbatively. We will also discuss the perturbative renormalization of the moving- 
NRQCD action from high-/3 Monte Carlo simulations. 
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Figure 1: Kinematics of B ^ ;r/v semileptonic decay, ps is the momentum of the B meson and p-,i is the 
momentum of the pion, q is the momentum transfer. 



Precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix is crucial for high 
precision tests of the standard model. The matrix element \Vub\ can be determined from studies 
of B ^ Tilv semileptonic decays in B-factories. Such decays, however, involve strong interaction 
dynamics also in the form of hadronic matrix elements; therefore precise determination of \Vub\ 
requires accurate calculation of these form factors. 

Lattice QCD provides a first principles non-perturbative approach for calculating the form 
factors f^{q^), fo{q^) in B nlv decays. Fig. |l] shows the kinematics of the process. Lattice 
simulations are problematic, however, when the pion has very large recoil momentum, i.e., when 
q-^ is small. This is because discretization errors increase as the pion momentum increases; and 
statistical errors become worse when the hadrons have large momenta. Lattice results are therefore 
available at the large q^ region only These simulations were done in the B meson rest frame, 
and two different techniques were used to handle the b quark on the lattice: non-relativistic QCD 
(NRQCD) and the Fermilab formalism. Because of the lack of data at small q^ error in \Vub\ is 
dominated by error in lattice results. 

This problem can be solved in two stages. First instead of working in the B meson rest frame 
we choose a moving frame in which the B meson is moving in the opposite direction of the pion. 
This reduces the momentum of the pion, but at the expense of increasing finite-lattice-spacing 
errors for the now high-momentum B meson. The increased errors in the B system are mostly 
associated with the b quark since most of its momentum is carried by the heavy quark. We work 
around this problem by using the moving- NRQCD (mNRQCD) formalism [|, ^ for the b quark. 

Moving-NRQCD provides an accurate discretization of a heavy quark moving at arbitrary 
velocity v. It is useful to parameterize the b quark momentum as pt, = unih + kwheie u = ps/rriB = 
{Es/tnBjv) is the 4-velocity of the B meson. The residual momentum k is much smaller than pB 
and ~ 0{Aq^jj). In mNRQCD one treats the velocity v exactly and discretizes the residual 
momentum only. The leading discretization errors therefore become 0{a^k^) instead of 0{a^p\), 
which could be as large as several GeV^. We refer the readers to Refs. ^ for details of the 
formulation. 

The mNRQCD action can be derived using the FWT transformation or by boosting the NRQCD 
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action to a moving frame [Q, ^. In this project we use an action which is complete through 0(1 /M) 
and tadpole improved (the tadpole factor mq is implicit in the equations): 



Evolution equation 

2 J \ 2n J ''''y 2n J y 2 J 



G^v+i= 1-^ 1-^ uIj l-^] 1-^ Gv. 



Action 



2Mo7 (1.1) 
5H= ^—a-B^+(-^ + -] (vA^)'-— ^(vA±)a(2) 

+ ^ y v,AtA±Ar + -J— y Af ) - (v • A±) Vv.A+Af Ar 



where 7=1/ Vl — v^, Mq is the input bare mass, n is the stability parameter, and B,„ is the magnetic 
field in the moving frame (we use 0{a^) improved electric field E and magnetic field B). The 
corresponding dispersion relation is given by 



E{k)=Eo + ^ Pl,+Ml^=Eo + ^ {ZpPo + kf + ZlMl (1.2) 

where Pq = yM^v. Here £0 is the zero energy, is the mass renormalization and Zp is the external 
momentum renormalization, which is a new feature in mNRQCD. These renormalization constants 
have perturbative expansions 

Eo{v) = ci^Ea{v)a, + (1.3) 

Zm(v) = 1+Ci,zm(v)a, + ..., (1.4) 

Zp(v) = l+ci,zp(v)a, + ... (1.5) 

Notice the dependence on the velocity. In the next section we will discuss how to exti"act the 
perturbative coefficients c\,C2 ■ ■■ from high-j3 Monte Carlo simulations. 

2. PERTURBATIVE EXPANSIONS FROM HIGH-/3 (WEAK COUPLING) 
SIMULATIONS 

Knowing how the bare action parameters renormalize on the lattice is crucial for accurate 
simulations. The first order coefficients (quenched with the Wilson plaquette action, 5H = 0) have 
been computed in standard lattice perturbation theory using diagrammatic techniques [^. Analytic 
calculations, however, are very challenging because of the proliferation of diagrams due to the 
lattice cutoff. This is particularly true at higher orders. An alternative is to use numerical methods. 
Here we determine the first order coefficients from high-j8 simulations [Q] and compare the results 
to the analytic calculations. 
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Figure 2: (a) Zero energy £'o(v) as a function of The volume is 6^ x 12. (b) 1/L infinite volume 
extrapolation of the first order coefficients ci £o(v). 



We start by doing ten sets of simulations, at volumes V = L^ xT = 6^ x 12, (8^, 10^, 12^) x 16, 
with j8 -values range from 10.0 to 60.0. The corresponding are of the order of 0.01-0.1. It is 
important to work at small couplings so that the theory enters the perturbative phase. Twisted 
boundary conditions are also used in the x and y directions to suppress non-perturbative effects due 
to Z(3) phases [^. In each simulation we measure the heavy quark energy at several momenta 



and fit the results to Eq. |L2| to extract the renormalization constants. We were able to work at 
small momentum, < 0.3, by using generalized boundary conditions for the quark field in the 
z direction Finally to obtain the perturbative coefficients we fit the renormalization constants 
obtained at different j8 to a polynomial in a^. 

Fig. ^ plots the zero energy £'o(v) against at several velocities at volume 6^ x 12. The 
"slopes" of the lines are the first order coefficients ci^eo{v) and the "curvatures" are the second 
order coefficients. Fig. |2|b shows the 1 /L infinite volume extrapolation of the first order coefficients 
ci,£o(v); and the final results are given in Fig. ^ Our numbers agree very well with the analytic 
results (given by the curve in the graph). We have also checked that reversing the fit order, i.e., 
extrapolate to infinite volume at each j3 and then fit to an expansion in a^, gives consistent results. 

Fig. ^ shows the external momentum renormalization Zp{v) obtained at 6^ x 12 and 12^ x 16. 
Despite the large statistical errors all the data are consistent with 1. We can conclude that there is 
little or no renormalization of the bare external momentum. This is consistent with the fact that 
V is protected from renormalization because of approximate re-parameterization invariance of the 
action 



3. B MESON DECAYS ON MILC FINE LATTICES 



In this section we present results for B meson correlators on MILC fine lattices. We study the 
dependence of B and B* masses on v. We also show how to extract the momentum renormalization 
from simulations. 
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Figure 3: Zero energy first order perturbative coefficients ci £o(v) after infinite volume extrapolation. The 
curve shows the analytic results obtained in Ref . [0] . 
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Figure 4: External momentum renormaUzation Zp(v) as a function of Us- Results are shown for volumes 
6^ X 12 and 12^ x 16. 



Simulations were done with the 0{\/M) accurate tadpole improved (mq = 0.8541) mNRQCD 
action at three different velocities v = (v,:,0,0), Vv = 0.0, 0.25 and 0.5. The MILC fine lattices 
have a lattice spacing of 2.258GeV-\ which is set by the T 2S-1S mass splitting [Q] ; the sea quark 
mass is uoamf^ ^ = 0.0062,0.031. The bare mass is set to aMo = 1.95 {n = 2) for the heavy quark; 
and we use the asqtad formalism for valence light quarks with a quark mass aniq = 0.031. It has 
been shown that, from previous NRQCD studies on heavy-light systems and the T spectrum, these 
quark masses reproduce the correct S., meson mass at v = 0. Gaussian smearings with a radius of 
2a were applied to the heavy quark at both the source and the sink. Once again we measure the 
meson energy at different momenta, and fit the results to Eq. 1.2 to extract the zero energy Eq{v), 
kinetic mass mkin{v), and momentum renormaUzation Zp{v). Thirty nine different momenta, with 
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Figure 5: (a) Bs meson mass mi^i„ as a function of bare velocity, (b) B^-B* energy splitting as a function of 
V. The experimental values are given by the horizontal lines. 



0<k^ < 0.45, were used. 

Fig. |5]a plots the mass as a function of bare velocity v^. Result at v = agrees well with 
the experimental value as suggested in the last paragraph. Our data show almost no dependence on 
the velocity, although a small increase can be observed as Vx becomes large. This suggests that one 
may need to tune (reduce) the input bare mass Mq as v increases. 

The energy splitting of Bg-B* is plotted in Fig. ^5. Again results show only a small velocity 
dependence. Data at v = 0, however, is about one a lower than the experimental value. It is 
believed that radiative correction of the a • B term in the action is responsible for this deviation 
from the experimental data. 

Finally the external momentum renormalization Zp(v) for Bg meson is given in Fig. ^ All 
the data points are consistent with 1 . This agrees with the results obtained in perturbation theory 
(see Fig. Qin the last section), and again reflects the re-parameterization invariance property of the 
action. 

Fig. m and ^ both reveal that statistical errors increase significantly between Vx = 0.25 and 0.5. 
It is therefore advisable to use a more advanced smearing technique, such as the so-called "magic" 
smearing [||], in future simulations. 



4. CONCLUSION 



Moving-NRQCD provides a promising approach for studying semileptonic decays at high 
recoil on the lattice. This project represents a major step forward towards this goal. We have 
demonstrated that perturbative expansions can be computed in high-j3 simulations with much less 
effort. First order results agree very well with the results obtained using analytic methods. Work is 
currently underway to extract higher-order coefficients from the data. In the second part of the paper 
we present results for B meson correlators on MILC fine lattices. A variety of quantities have been 
measured at different B meson velocities; and we observe little or no velocity dependence. This 
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Figure 6: Momentum renormalization Zp{v) of By as a function of input velocity. 

is significant since it indicates that we do not need to tune the input bare parameters as v changes. 
On the other hand our results show a sharp increase in statistical errors as v becomes large. This 
suggests that it may be difficult to reach a very low value. It is, however, not necessary to give 
the form factors at very small since the main experimental focus is on 8GeV2 <q^< 16GeV2. 
Simulations for B nlv decays will start soon and we expect to have new results in the next year. 
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